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Abstract 

The Regge-Wheeler equation for black-hole gravitational waves is analyzed for large negative 
imaginary frequencies, leading to a calculation of the cut strength for waves outgoing to infinity. 
In the — limited — region of overlap, the results agree well with numerical findings [Class. Quantum 
Grav. 20, L217 (2003)]. Requiring these waves to be outgoing into the horizon as well subsequently 
yields an analytic formula for the highly damped Schwarzschild quasinormal modes, including the 
leading correction. Just as in the WKB quantization of, e.g., the harmonic oscillator, solutions in 
different regions of space have to be joined through a connection formula, valid near the boundary 
between them where WKB breaks down. For the oscillator, this boundary is given by the classical 
turning points; fascinatingly, the connection here involves an expansion around the black-hole 
singularity r = 0. 
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I. INTRODUCTION 



Black-hole axial gravitational waves of angular momentum £ > 2 are described (in units 
c = G = 2M = 1) by the Regge- Wheeler equation (RWE) [1, 2] 
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where x = r + ln(r— 1) is the tortoise coordinate and r the circumferential radius; V(r) 
accounts for the Schwarzschild background. The long-range nature of this potential, 
V(x) — £{£+!) jx 1 ~ 2£(£+l) lnx/x 3 for x — > oo [3], is well known to cause a branch cut in 
the (retarded) Green's function of (1) on the negative imaginary axis (NIA) in the cu-plane. 

For 7 = iuj I 0, this cut causes an algebraically decaying late-time tail in the gravitational- 
wave signal [4, 5]. For moderate 7 and 2 < £ < 4, it has recently been investigated numeri- 
cally [6], leading to a clear conjecture for the large-7 behavior. In a separate development, 
the strings of quasinormal modes (QNMs) parallel on both sides of the cut and close to it 
seem to offer clues to the quantum theory, in particular to a calculation of the Bekenstein 
entropy in loop quantum gravity and to the quantum of area [7]. These motivate studying 
also the branch cut asymptotically, which will turn out to have ample independent interest. 

In each ^-sector and in the frequency domain, the RWE is one-dimensional, so the above- 
mentioned Green's function G(x, y) = G(y, x) can be written as 

G(x, y; u) = — y < x . (3) 

J{0J) 

Here, / solves (1) with the left outgoing-wave boundary condition (OWC) f(x— 00, uo) ~ 
1 • e~ tuJX and thus represents waves going into the horizon, while g(x^oo, iv) ~ 1 • e tuJX corre- 
sponds to waves going to infinity; J = gf — fg' is their Wronskian. In the physical region 
iGK, these asymptotic definitions are unambiguous only for Imu > 0, from where the func- 
tions are continued analytically. The normalizations of /, g have been fixed for definiteness, 
but J in the denominator renders (3) normalization-independent. Since V(x^— 00) ~ e x , 
the function f(x,cu) is single-valued in uj. Thus it is intuitively clear, and readily shown [6], 
that the branch cut in G can be expressed in terms of the one in g. 

Focusing on the latter, we define g±(w) as the continuations from Re±w > 0, and Ag = 
g + — <?_. Since g±(x, —27) ~ 1 • e lx satisfy the same (linear, second-order) wave equation, 
Ag(x, — ry) ~ Q-e' x is the small solution oc g(x, The simple symmetry g{— uj*) = g*{u) 

renders Ag imaginary, so we introduce the real cut strength q through [4, 8] 

Ag(x, -ij) = iq(i)g(x, . (4) 

Since g is defined by the OWC at x — > 00, (4) defines q not merely x-independent, but 
rather independent of V(x) at any finite x: if, say, V\{x>L) = V2(x>L), the corresponding 
qi and q 2 are identical. Thus, q economically characterizes Ag (and ultimately AG), and 
our task can now be specified as calculating the asymptotics of ^(7^00). 

II. WKB SOLUTIONS 

For Imw < 0, in particular on the NIA, the simple asymptotic definition of g(ou) becomes 
all but meaningless, since it is impossible to distinguish the decaying component 
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(to be set to zero) from the 'pure' outgoing wave ~ e lx , where this limiting form however 
has algebraic corrections to all orders in V . One way out is to complement the analytic 
continuation in frequency with one in space [8-10], so that the product iux retains a negative 
real part. However, in terms of x, the very equation (1) is multiple valued, so that the 
analysis henceforth will proceed in the complex r-plane, viz., 



[r 2 (r-l)X + r(r-l)d r - (r-l){£(£+l)r - 3} + u 2 r A ]g = . 



(5) 



It is possible to impose the OWC for g+(r, — ij) [g~(r, — ry)] directly and stably for 
r — > — oo, and continue the solutions to the physical r > 1 in the upper (lower) half 
plane [8]. That is, apart from a trivial overall phase, g± are the same solutions as r grows 
from — oo, until they are prescribed to encircle the singularity r = in opposite directions. 
Hence, closer study of this point should shed light on their difference Ag. 

At least away from the singularities r = 0, 1 and the anti-Stokes lines (v.i.), one expects 
to have asymptotic expansions g a (r,u>) and g a (r, -co), with 



gjr,u) ~ [(r-l)e r f 
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where the first factor is just a plane wave in the tortoise coordinate. Substitution yields 
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where N = 8(^ 2 ) = with v = ^(£—l)(£+2). That is, these two orders still agree 

with the upshot of cu-expanding the standard WKB expression, though this is no longer true 
for <? 3 , which involves V (of course, there are higher-order corrections to WKB as well). 

The anti-Stokes lines of the RWE on the NIA for this expansion are shown in Fig. 1. They 
are the curves where |(r— l)e r | = 1 (i.e., Rex = 0); for 7 — ■> 00, these are the boundaries 
where the solutions (6) change between exponentially growing and decaying characters. 

Thus, the solution which is small for r — > — 00 can be continued to the region including 
x > 0, but this only yields the exponentially growing part one knew all along. From there, the 
solution cannot be continued back to x < (where its imaginary part could be identified) 
because of the Stokes phenomenon. For the latter continuation, we apparently have to 
pass through the eye of the storm — the black-hole singularity r = [11]. Near r = 0, 
the expansion (6) is not valid, for higher-order terms only become small if |r| ^> 1/^/7. 
One therefore has to match (possibly different) expansions in terms of g a in the regions 
|argr — 7r| < j and |argr| < | respectively. This is similar to the connection procedure 
near classical turning points in Bohr-Sommerfeld quantization. However, near the latter 
it is merely the asymptotic expansion which breaks down, while the original (Schrodinger) 
wave equation is perfectly regular. In our case, the connection has to be carried out across 
a singularity of the RWE. 
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III. CONNECTION FORMULA 



Series expansion around r = is standard [12]: 

^(r) = r 3 + 6 ~y +1) r 4 + G(r 5 ) , (9a) 
5 

, / x i 2v N £(£+l)N 2 fi 2 - , ^, 4 , 

^ 2 ( r ) = r- 1 + _ + — r + r 2 + ^(r) lnr + C(r 4 ) , (9b) 

where both error terms are single-valued, and where we introduced the algebraically special 
frequency Q = —iN/2 [8, 13]. However, because of the large term oc u 2 in the RWE (5), 
higher-order terms in these expansions can only be omitted if \r 2 uo\ <C 1 (cf. the given terms 
of ip2)- Thus, there is no overlap with the region of validity of g a , where matching could be 
carried out. Still, the above series for ipi^ will be useful for comparison, inter alia yielding 
the exact monodromy [14] ip 2 {re 2m ) = fair) + i^(Q 2 — uj 2 ) ipi(r) and ifj 1 (re 2m ) = ipi{r). 

A usable matching solution follows by effectively resumming the large-7 parts of higher- 
order terms in (9). In practice, it is more convenient to set r = t/^y and sort powers 
of 7: 

[t 2 d 2 -td t -3- t 4 ]^ = — [2t 3 d 2 - t 2 d t - {£{£+!) + 3}t]^ + 0(7^) = R(t) , (10) 

^ = ^(°)+^( 1 ) + ... . (11) 
The lowest order follows by equating the rhs of (10) to zero: 

^« = ^W|). (12b) 
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These have not been written in terms of modified Bessel functions, since the subsequent 
matching is best done on the anti-Stokes lines where t 2 /2i is real. Subsequently, ipf^ figures 
as an inhomogeneous term in the equation for solved by 
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Note that the occurrence of Y 1 does not spoil the analyticity of ip^ . In particular, for t — > 
the above readily reproduces the t 4 -term found in (9a) by direct expansion. The counterpart 
for ip2 reads 
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where we introduced a 'counterterm' in order to keep the second integral finite near s = 0. 
Again, one verifies that the leading small-t correction [oc t° in (9b)] is reproduced correctly. 
Using the standard branching properties of the Y n one now finds that, up to this second order, 
ip 2 (te 2m ) = ^(O + if l 2 i>i{t)-, so that the exact monodromy is approached for \uo\ 3> 

In fact, we can follow the transformation of tp^ under rotation in more detail by using 
the standard J n (-z) = (-) n J n (z) and Y n (e ±i7T z) = (-) n [Y n (z) ±2iJ n {z)}, implying 



^(it) = -i^\t), 
•Xo) ( . A _ • ,,(o), 



(17a) 
(17b) 

(17c) 
(17d) 



These relations streamline the asymptotic expansion. Namely, look along argt = |, where 
ipfl e" 37 ™/ 4 G R and ip^l G R, so that their expansions will read 
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The occurrence of (9(t 3 / v /7) [which however does not generate additional t°/ Ai /7-terms] means 
that these hold for I <^ t 7 1 / 6 or 7 -1 / 2 < r < 7 -1 ^ 3 , conveniently handled as a double 
asymptotic expansion in t and A = 7/t 6 . The power-law asymptotic corrections in t follow 
directly from the RWE; the real nonlocal information is contained in a, [3. 
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We can use the general rules (17) to transform (18) to argt = and demand consistency 
for the (oc e~* / 2 ) part that dominates for | < argt < This leads to a G R and 
Im p = —a. By subsequently demanding consistency also for the combination of ip\ and 1^2 
which is minimal in the same sector, one finds 

= -(2+i)a . (19) 
Thus, this algebraic exercise circumvented directly expanding the integrals (15) for 



IV. MATCHING 



For |argt — tt\ < one can do the analogous expansion of g+(t, —27) = g a (t, —27) in (6). 



Again, one finds a Gaussian form if t <C 7 1 / 6 [15]: 

3 15 t 3 
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the first factor comes from encircling the horizon r — 1. The leading corrections verify the 
consistency of expansions (18) and (20), obtained in very different ways; all that matters for 
the matching is [•••] = 1 + • £%/7 + • • • ■ 
Comparison shows that 

g + e-**> = ^(3* 7 3/2 + (2-3*)a 7 + 0(Vt))^i + + ^ + 0(7 _S/2 ))lfe , (21) 

which can be matched back to solutions in terms of g a on argt = |, yielding 



g+i-i-y) = g^-i-y) + 2e 2 ^ 



(22) 



in the region (bounded by anti-Stokes lines) including r — \. Hence, in particular in the 
physical part x < of that region, one has 



Ag(— 27) = 2i Im <?+(— 27) = Ai 
= iq{l)g{+il) . 



a 
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(23) 



Already, the numerically observed asymptotics [6] of 9(7) have been confirmed. For the 
corrections, it remains to calculate a in closed form. Let us start with the //-dependent term 
in (14), which clearly cannot be combined with the other two. Straightforward manipulations 
yield its contribution to ip^ as (argt = j) 



-v 
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By considering the asymptotics of ^(dz/y/z) Ji(z)h[ 1 \z) (H are Hankel functions) in 
the upper-half A"-plane [again using the formula for Y\{e m z) above (17)], one convinces 
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oneself that in fact f£°(dz/y/z) Ji(z)Y 1 (z) = — f^°(dz/\/z)J 2 (z), which is also necessary for 
this contribution to a to be real. The latter integral is tabulated as J^°(dz/y/z) J\ [z) = 
r(i) 4 /127r 5 / 2 [16]. Thus, the present contribution to a reads 

The last term of (14) analogously leads to the integrals / °° dz y/z J\{z) Jq(z) = 
— J Q ° C dz yfzY\{z) Jo(z) = r(l) /167T 5 / 2 , for a contribution 

rfi) 4 



However, the first term of (14), with its higher power of s, leads to diverging integrals: 

" K dzz 3 / 2 J*(z) ~ ^K 3 ' 2 + ^-VKcos(2K) + c 3 + 0{K~ x l 2 ) , (27) 

j 07T 2lX 

r K 1 

/ ^2 3 / 2 Jx (z)Y 1 (z) y/Ksm{2K)+d 3 + 0{K- x l 2 ) . (28) 

Jo 2tt 

^/ ^3 — — c 3; these lead to a real contribution 

a 3 = -27TC 3 . (29) 

Unfortunately, the general f*dzz 3 / 2 J 2 (z) = (X 9 / 2 /18) 2 F 3 (|, f ; 2, 3, f ; -fT 2 ) does not help 
directly, since not enough seems to be known about the asymptotics of 2^3- Instead, one 
can proceed as follows: d 3 = — c 3 can again be proven by studying J* K dz z 3 ^ 2 Ji(z) h[ 1 \z) in 

the upper-half plane. It is then logical to also consider J A dz z 3 ^ 2 H[ 1 \z) 2 , in which one can 
take K — > ioo. One finds 

* = ^f W/2 *>M = SS' (30) 

The rest is straightforward: (25), (26), and (29) with (30) can be added and substituted 
into (23), from which one can read off our final answer 

9(7) ~ 4 cos(27r 7 ) + ^ [1 - 1(1+1)] + O^ 1 ) , (31) 

where rQ) 4 /127r 3 / 2 = 2.586... The result for the leading term has appeared before in 
Ref. [6]; preliminary results from a transmission- amplitude calculation seem to support the 
form for the corrections, in particular the ^-dependence [17]. 



V. DISCUSSION 

The above analysis, for r, 7 — > 00, is the third instance where the mathematics of the 
RWE near r = has been seen to influence the goings-on in our universe r > 1. The first 
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instance is the question of 'anomalous' vs 'miraculous' waves outgoing into the horizon, i.e., 
r | 1 and 27 e N ([8], Section VI). The second instance is the dynamics at and around the 
algebraically special frequency 7 = iV/2, where the RWE has closed-form solutions, whose 
global behavior in the r-plane can therefore be traced [6, 8]. 

Besides the axial waves described by the RWE (l)-(2) there are also polar waves, de- 
scribed by the Zerilli equation [18]. However, the latter's 'intertwining' or supersymmetry 
relation to the RWE yields its cut strength as £(7) = [(N/2--f)/(N/2+-f)]q(-f) [8]. Thus, q 
and q agree up to an overall sign, plus Oij^ 1 ) corrections immaterial to (31). 

For an outlook, the first obvious item is the numerical verification of (31), in particular 
of the sub-leading correction. A brute-force increase of numerical precision in the existing 
method is unlikely to suffice, especially for i > 3. More promising is to calculate g±{—i r y) 
directly on the NIA, instead of through extrapolation to this axis. In the series for g(u>) 
[6, 12], the problem for Re ui — > is not convergence but rather obtaining the individual 
(irregular-hypergeometric) terms reliably. Working this out should be mainly a matter of 
time, but it remains to be seen if it would sufficiently extend the range of validity in 7. 

Related to this, it is worth re-emphasizing [8] that the present method involves two 
conceptually separate steps: (a) the continuation in r, stabilizing the OWC at infinity 
in the lower-half a;-plane, leading to a well-defined computational problem, and (b) the 
asymptotics, by which one can actually solve this problem analytically for large 7. The 
second step is optional [9], and numerical integration of g + from r = —00 (combined with 
standard series solutions for /) should soon open up the third u;-quadrant (i.e., behind the 
cut) to direct exploration, especially for at most moderate damping. 

It would also be interesting if this work could be compared to the closed-form ex- 
pression for 5(7) in (31)-(33) of Ref. [4], through the coefficients d^} of an expansion 

9 SlL-oo d^UL+v', here, ul+ v are Coulomb wave functions. As usual, the d^ satisfy 
a three-term recursion relation (given in a simplified form, which can be made purely real 
on the NIA, in Section VI F of [12]); < v < 1 is to be chosen such that d^ is the min- 
imal solution to this relation for both L — > — 00 and 00. The adiabatic Ansatz of slowly 
varying d^^/d^ readily yields asymptotic solutions for 7 — > 00, except near L = ±7 and 
L = 0. Following Ref. [19], one can try to develop connection formulas near these three 
points, which would determine v analytically. At the 'turning points' L = ±7, the three- 
term recursion asymptotically degenerates into a two-term one, and the connection proceeds 
exactly as in [19]. Near L = 0, however, all three terms are of the same order (in 7), and the 
recursion remains in the form of analytically intractable continued fractions (i.e., without 
simplifying to products). Thus, this route for now seems unfeasible. 

At least to leading order, the method of continuation through the vicinity of r = has 
meanwhile independently been used to calculate the high-damping QNM frequencies [11]. 
On the one hand, this yields more information on the QNM wave functions than the 
continued-fraction technique [19]. On the other, this clearly establishes a relation to the 
present problem of the branch-cut strength 9(7). Thus one can expect also numerical re- 
sults for these two high-damping aspects of the RWE to bear on each other, and an extension 
to other black-hole models. Conversely, this work's progress on the corrections should help 
finding the (9(|c<j| _1 / 2 ) terms for the asymptotic QNM frequencies. Meanwhile, this has 
indeed been possible; see the Appendix. 

In particular, these developments render it urgent to study the branch cut of the RWE's 
Green's function also for other values of the spin s, where (2) corresponds to s = 2. Notably, 
the highly damped electromagnetic (s = 1) QNMs are predicted [11, 19] to approach the 
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NIA asymptotically so that one expects an even closer, though yet unknown, relation to the 
branch cut in that case. 
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APPENDIX A: HIGH-DAMPING QNMS 



For frequencies close to the NIA, it should certainly be possible to analytically con- 
tinue (22) and find those 7 for which g + also satisfies the OWC into the horizon, i.e. 
g+(ou) oc f(cu). These should identify the highly damped QNMs. Here, / is characterized by 
its monodromy around the horizon: 



f((r-l)e? m ,u) = e?™f(r-l,u). 



(Al) 



Inspecting (6), it seems that 1 fiyj) = g a (~ w)\ but this is deceptive. While it is true that 
f{—oj) = g&(uj) in the lower-half w-plane, the asymptotic nature of the large (near r = 1) 
solution g a (—ou) means that naive rotation of r — 1 only confirms the desired monodromy to 
dominant order, whereas / is required to obey it exactly. 

The solution is to, in the spirit of [11, 17], do the rotation along the anti-Stokes contour 
in Fig. 1, where neither solution dominates the other (skipping r = on the inside). By the 
normalization of / one knows the dominant component, so we take the Ansatz 

f(-il) = gSl) + c(7)0aH7) • (A2) 
Continuing this from the physical x = to <C |r| <C 1 with argr = |, one finds 



/(-i 7 ) ~ e~ l7r V 2/2 + c( 7 )e l7r7 e~* 2/2 
" (i-2)a 
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(A3) 



the second line followed by comparison with (18). As before, (17) makes quick work of con- 
tinuing this to argr = — |, where it can be matched back to a combination of ^(±^7). One 
thus finds f((r-l)e 27T \-ij) = e~ 2ni ~< g^r-1, i 7 ) + [2i(l-a/^) + 0(7) e 2n ^]g a (r-l, -i 7 ). In- 
deed / as in (A2) obeys (Al) to dominant order for any c, while the subdominant (oc c) term 
by itself has the opposite monodromy, corresponding to incoming waves — both as stipulated 
above. Equation (Al) holds exactly for 



c( 7 ) 



1 - «A/7 
sin(27T7) 



(A4) 
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As required, c( 7 GM) G M., since the exponential tail of V(x^— oo) does not generate a 
branch cut in f(ou). Also the poles for 27 G N are not surprising, since the RWE is known 
to have such anomalous points for 27 = 1, 2, . . ., with exactly one exception at 27 = N [8]; 
clearly, the latter is beyond the reach of the present asymptotics. Expressing f(±u), g(±cu) 
not through g a but in terms of each other, the S-matrix could be read off; cf. [17] for the 
leading order. 

Combining (22), (A2), and (A4), all that remains is to asymptotically solve 
= J(-i 7 ) oc sin(27r 7 ) + 2e 2 ™ 7 (i + (\ - 

« sin(27r 7 ) + 2*e 2 ™ 7 (l - • (A5) 

Reexpressing the answer in terms of uj and substituting a from the main text, one obtains 



where the prefactor of the correction evaluates to 0.097007. . . The u n have automatically 
come out on the physical sheet of g + ; contrast Ref. [19], which is not sheet-specific. Agree- 
ment with (29)-(30) in Ref. [20] is excellent; as anticipated in the Discussion, this numerical 
confirmation of the value of a via the QNMs greatly supports (31) for the cut. To my 
knowledge, this is the first time that the correction to the highly damped Schwarzschild 
QNM frequencies has been calculated analytically. 
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